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Ground-state properties of gapped graphene using the random phase approximation 
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We study the effect of band gap on the ground-state properties of Dirac electrons in a doped 
graphene within the random phase approximation at zero temperature. Band gap dependence 
of the exchange, correlation and ground-state energies and the compressibility are calculated. In 



y—' , addition, we show that the conductance in the gapped graphene is smaller than gapless one. We 



also calculate the band gap dependence of charge compressibility and it decreases with increasing 
the band gap values. 

d . PACS numbers: 73.63.-b, 72.10.-d, 71.10.-w, 73.50.Fq 



I. INTRODUCTION 



j^ ' Graphene is a flat monolayer of carbon atoms tightly packed into a two-dimensional (2D) honeycomb lattice and it is 

^ ' . . . 

I , a basic building block for all nanostructured carbon. This stable structure has attracted considerable attention because 

f~^ of experimental progress''^ and because of exotic chiral feature in its electronic properties and promising applications^. 



Very recent experiments on both a suspended graphene and a graphene on substrate have found remarkably high 
mobility 2 x 10^ cm^/Vs for carrier transport^ at room temperature which is two order of magnitude higher than the 



1-^ , mobility of silicon wafer used in microprocessors!^ 

Q ' An interesting feature of graphene which makes it very applicable in semiconductor technology is the opening a 

gap in the band energy structure of graphene. There are several scenarios to open a gap in the band energy structure 



CO of graphene. One is the finite size effect by using graphene nanoribbons, where gaps give rise to constriction of the 
OO , electrons in the ribbon and it depends on the detailed structure of ribbon edgesi^ More precisely, the band gaps with 



armchair shaped edges originate from quantum confinement and the value of the gap depends on the width of the 
ribbon. For zigzag shaped edges, on the other hand, the band gaps arise from a staggered sublattice potential due to 
magnetization at the edgesi^ Importantly, edge effects plays a crucial role in transport properties. The gap engineering 
can be also achieved through doping the graphene with chemical speciesZ due to the translational symmetry breaking. 
The electronic properties of a graphene interacting with CrOa molecules has been calculated by using ah initio 
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H calculations.?' This type of calculations, predicts opening a gap about 0.12 eV at the Dirac point. Another scenario 
is graphene by placing it on top of an appropriate substrate which breaks the graphene sublattice symmetry and, 
therefore generates an intrinsic Dirac mass for the charge carriers.^ Typical substrate is made of hexagonal SiC with 
a gap about 0.26 eV. A recent band structure calculation for a graphene on top of a hexagonal boron nitride crystal^^ 
has been shown a band gap about 53 meV. The gap can also be generated dynamically by applying a magnetic 
fieldjii Moreover, when both mono-and bilayer graphene material are covered with water and ammonia molecules, a 
gap induce in the spectrum of energy.^^ Interestingly, the mechanism that electrons hopping on a honeycomb lattice 
with textured tight-binding hopping amplitudes, the Kekule texture, generates a Dirac gap J^ Eventually, It has been 



suggested that a small gap can be opened on the Dirac points due to spin-orbit coupling or Rashba efFect^^ which 
makes the system a spin Hall insulator with quantized spin Hall conductances^^ 

Recently, the local compressibility of graphene has been measured^^^ using a scannable single electron transistor. The 
measured compressibility is claimed to be well described by the kinetic energy contribution and it is also suggested that 
exchange and correlation effects have canceling contributions. From the theoretical point of view, the compressibility 
was first calculated by Peres et al^ considering the exchange contribution to a noninteracting doped or undoped 
graphene flake. A related quantity dfi/dn (where /i is the chemical potential and n is the electron density) is recently 
considered by Hwang et al}^ within the same approximation and they stated that correlations and disorder effects 
would introduce only small corrections. This statement is only true in quite large density doped values. To go 
beyond the exchange contribution, the correlation effects were taken into account by Barlas et alr^ based on an 
evaluation of graphene's exchange and random phase approximation (RPA) correlation energies. Moreover, Sheehy 
and Schmalian^S by exploiting the proximity to relativistic electron quantum critical point, derived explicit expressions 
for the temperature and density dependence of the compressibility properties of graphene. Importantly, the effect 
of disorder and many-body interactions on the compressibility has been recently studied by us.~ We successfully 
demonstrated the importance of including correlation effects together with disorder effects in the thermodynamic 
quantities. It should be noticed that all these theoretical efforts have been carried out for a gapless graphene. 

Our aim in this work is to study the ground-state properties in the presence of Dirac gap and electron-electron 
interactions. For this purpose, we derive the gap dependence of the dynamic polarization function for a doped 
graphene to calculate the scattering rate, ground-state energies and the compressibility of the system at the level of 
RPA including the opening gap at Dirac point. 

We consider different on-site energies for atoms in two sublattices in graphene which is established experimentally 
to be important when an appropriate substrate such a boron nitride or SiC is used. The compressibility decreases by 
increasing the band gap values due to the sublattice symmetry breaking. 

The rest of this paper is organized as follows. In Sec. II, we introduce the models for dynamic polarization function 
and ground-state energy calculations . We then outline the calculation of d.c conductivity and compressibility. Section 
III contains our numerical calculations of ground state properties. We conclude in Sec. IV with a brief summary. 

II. THEORETICAL APPROACH 

We consider a Dirac-like electron in a continuum model interacting via a Coulomb potential e^/e?' and its Fourier 
transform Vq = 2Tre'^/{eq) where e is the average background dielectric constant (for instance, e ~ 5.5 for graphene 
placed on SiC with the other side being exposed to air) having an isotropic band gap at Dirac points. If one assumes 
that the sublattice symmetry is broken and Ua, ctb are on-site energies of atoms A and B, respectively, then the 
contribution of the on-site energies in Hamiltonian of graphene can be writter>ii as 



Hi = ^[aaat(k)a.(k)+a66^(k)6,(k)] 

k,cr 

= 5]#^[a+r"-a_r3]®a3vPk,. (1) 

k,<T 

where a+ = (a^ + at,)/2 that corresponds to the same carrier density on two sublattices, a_ = («„ ^ ctb)/2 is carrier 
imbalance on two sublattices that leads to break the inversion symmetry and r° is 2 x 2 unit matrix, r^ is a Pauli 
matrix that acts on Kj^- and K^ two-degenerate valleys at which tt and tt* bands touch and a^ is Pauli matrix that act 
on graphene's pseudospin degrees of freedom. Consequently, the nonintcracting Hamiltonian for a gapped graphene 
is given by ^o = Ek.<T *k o-^o*k,^ where 
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where k — k^ + iky and a is the spin of charge carrier. Here, v — 3ta/2 is the Fermi velocity, t is the tight- 
binding hopping integral, a is the spacing of the honeycomb lattice. For the hexagonal crystal structure of graphene, 
a = 1.42A is the carbon-carbon distance, the tight-bonding hopping energy is i = 2.8 eV and the bare Fermi 
velocity is w = 10^ m/s. In the nonintcracting Hamiltonian, Ho the reference energy {aa + af,)To/2 is subtracted 
and the energy gap is defined as 2 A = (a;, — aa) where we expect A < t. The corresponding four components 
pseudospinor of the nonintcracting Hamiltonian is "^l.^. — ( V-'+,o-, V'+,<t: '0-,o-j''/'-,<t )■ It is easy to diagonalize the 
nonintcracting Hamiltonian based on pseudospinors in the conduction and valance band of energies with eigenvalues 



given by ±\/^^w^fc^ -I- A^. Importantly, the low energy quasiparticle excitations in a gapless graphene are linearly 
dispersing and it is valid for energy less than 1 eV. Accordingly, the validity of the nonintcracting Hamiltonian given 



by Eq. ([2]) to explore graphene properties is to the case which ^/h'^v^k^ + A^ < leV. On the other hand, we shall 
achieve to a conventional two-dimensional electron gas system by setting hvk/ A <^ 1. 

Finally, the total Hamiltonian including the electron-electron repulsion interaction is given by 

where 5* is the sample area and TV is the total number operator. The presence of a neutralizing background of 
positive charge is explicit in Eq. ([3|). As we mentioned in introduction section, this kind of Hamiltonian can be 
used in graphene by placing it on top of an appropriate substrate that breaks the graphene sublattice symmetry and 
generates an intrinsic Dirac gap. 

A central quantity in the theoretical formulation of the many-body effects in Dirac fermions is the nonintcracting 
dynamical polarizability functioni^i^^i^^ x^^'' (Qi *^: M > ^) where /i is chemical potential. Here, we would like to 
emphasize that we have calculated the gap dependence of the nonintcracting polarization function for doped graphene 



however the vacuum polarization function in which /i = A has been calculated by Kotov et a?^. They studied 
the distribution of polarization charge induced by a Coulomb impurity for undoped graphene. However, we would 
like to study the ground-state properties for doped graphene sheets. To achieve this goal, we write the dynamical 
polarizability function in terms of one-body noninteracting Green's function 

(P'k f dw 



x(°)(q,f^,M) = -z 



-Tr[z7oG("' (k + q, c^ + f}, m)»7oG(°) (k, u, ^l)] 



(27r)2 j 27r 
where one-body noninteracting Green's function^^ by using the noninteracting Hamiltonian is given by 



(4) 



in which 7-matrices are related to Pauli matrices by a^ — ~ijo and a^ = {—l)^a^jj for j = 1, 2 and kp is the Fermi 
momentum related to the density of electron as given by kp = (47rn/g)^' ■^. g — g^ gs — 4: is valley and spin degeneracy 



and 9 is the Heaviside step function. The chemical potential is given by /i = yft^u^kp + A^ at zero temperature. 
After implementing G*-*^-* (k, uj, /i) in Eq. ([4]) and calculating the traces and integrals, the result is given by the follow 
expression 
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where Sq = Hvq. By setting A = 0, it is easy to determine that Eq. ^ reduces to the noninteracting dynamic 
polarization function of the gapless graphene sheetJ^ Furthermore, for the half-filed gapped graphene sheet, the 
noninteracting dynamic polarization function, vacuum polarization, is given by^ 
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Using the above results for the noninteracting polarization function on the imaginary frequency axis, the density 
of state at Fermi energy is calculated as 



D{eF) = D^ep) f(l + AV4)'^'1 0{ep) 



(8) 



where _D''(ei?) = g£-p/2TThv'^ is the density of states of gapless graphene?^ Note that we define e-p — hvk-p- The Hnear 
correction of expanded gapped polarization function for A is zero, however the quadratic correction is easy obtained 



X^ 'iq,iuj,fi) ~ x^ 'iq,i^,fJ') 



(0), 
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Where the exphcit expression of x^'^^ il: ^^i m) 



A=o 



is given by our group.-- Now, we are in the stage to use the 



nonintcracting polarization function given by Eq. ([6]) to calculate some physical quantities. 

A. Transport scattering time in a gapped graphene 

As a first application of the nonintcracting polarization function, we would like to calculate the gapped graphene 
transport scattering time by randomly distributed impurity centers in the relaxation time approximation.^ The valid- 
ity of the Born approximation is discussed by Novikov^^ and here we use this approximation to calculate qualitatively 
the graphene transport scattering time. To this purpose, the transport scattering time is given by Boltzmann theory, 
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where Vi{q) = ^^^^ exp(— gd) is the Coulomb scattering potential between an electron and an out of plane impurity, 
e{q) is the static RPA dielectric function appropriate for graphene, €{q) = 1 — VqX^^' {q,0, ^i), rii is the density of 
impurities and d is the setback distance from the graphene sheet and s^s' being ±. Since we consider large charge 
carrier density and elastic scattering, we can therefore neglect interband scattering process. F''(q, q + kp) is the 
overlap of states (/3 = ±), which can be easily calculated from the pseudospinors of Hamiltonian, Eq. ([2]). The result 
is as follow 
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where (f> is an angle between k and q. Graphene conductivity can then be calculated by the Boltzmann transport 
theory with a = {e'^ / h)2T{eF)vkp . The properties of graphene's Dirac fermions depends on the dimensionless coupling 
constant agr — ge^/veh. 

B. RPA ground state energy in a gapped graphene 



The ground-state energies is calculated by using the coupling constant integration technique, which has the contri- 
butions Etot = Ekin + -Ex + Ec- The kinetic energy per particle is easy calculated as 2£_f[(1 -f- A^/e|,)^/^ — A'^/el.J/S. 
The first-order, exchange contribution per particle is given by 
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To evaluate the correlation energy in the RPA, we follow a standard strategy for uniform continuum models — 

Since x^''^(q, iw, /i) is linearly proportional to q at large q and decrease only like u!~^ at large ui in both gapped and 



gapless graphene, accordingly the exchange and correlation energy built by Eqs. ^2]) and ^3]) are divergent . ^^'^^ In 
order to improve convergence, it is convenient at this point to add and subtract vacuum polarization, x*-'^^(q, iw, /i = 
A), inside the frequency integral and regularize the exchange and correlation energy. Therefore, these ultraviolet 
divergences can be cured calculating 

1 r J2 r+oc 

^e.^-— y^.J^ ^.^X^°)(q,..,.) (14) 

and 

where Sx'"^^ is the difference between the doped {ji > A) and undoped (/i — A) polarization functions. With 
this regularization, the q integrals have logarithmic ultraviolet divergences.— we can introduce an ultraviolet cutoff 
for the wave vector integrals kc = Akp which is the order of the inverse lattice spacing and A is dimensionless 
quantity. Once the ground-state is obtained the compressibility n can easily be calculated from n^^ — r? — g„2*°* j 
where the total ground-state energy per particle is given by Jstot = (^Ekin + ^Ex + fc|f^^. The compressibility of 
noninteracting gapless graphene is Kq — 2/{neF) and the compressibility of noninteracting gapped graphene is given 
by^o = 2/nef(l-KAV£|)i/2. 

III. NUMERICAL RESULTS 

In this section, we present our calculations for the ground-state properties of gapped graphene in a continuum model 
at low energy, using the model described in the previous section. Our results are considered a n-doped graphene sheet 
with agr — 1 and 2 where it is a typical value thought to apply to graphene sheets on the surface of a SiC or boron 
nitride substrates. We expect that the physics of gapped graphene is different from gapless graphene due to the 
sublattice symmetry breaking. In the follow, we will investigate these differences, quantitatively. 

Fig. 1 shows the noninteracting dynamic polarization function, x (q^ *'^j m) of both gapped and gapless graphene 
in units of the noninteracting density of state at the Fermi surface, D{eF) as functions of q/kp and lu/ef- In both 
cases, X {'ii *^i /^) linearly diverges with q at small wavelength region and decays as 1/uj at large frequency for finite 
A values due to interband fluctuations in contrast to the ordinary 2D electron gas. 

The function — x (Qi 0, fi) contains a number of noteworthy features is shown in Fig. 2(a). First, as we mentioned 
before, the q — > limit of the static polarization function is a measure of the number of excited states. Second, the 
derivative of x'^^Hli 0: /i) at q = 2kF is singular at finite A values the same as the normal 2D electron gas. Note that 



X*'°-'(l'Oi m) at a = is a smooth function. We stress here that the second order correction of the noninteracting 
polarization function is mostly responsible to this behavior. This singular behavior is responsible for several interesting 
phenomena such as Friedel oscillations and the associated RKKY interaction.^® Interestingly, x^"-* (l: Oi m) reduces to 
the contentional 2D noninteracting dynamic polarization function at very large A values. In this figure, we have shown 
^X (QjOj/i) at a = lOej^ which is exactly the same as the conventional 2D noninteracting dynamic polarization 
function up to mid q values. The behavior of x (kpiitjjM) in unit of the density of state of gapped graphene for 
various A is displayed in Fig. 2(b). 

As an application of the noninteracting polarization function given by Eq. ((6]), we calculate the electric conductivity 
using the Boltzmann equation, we assume d = lA and agr ~ 2. Band gap and density dependence of d.c. conductivity 
are shown in Fig. 3. Increasing disorder (increasing n^ or decreasing d for charge-disorder potential ) decrease the 
a however increasing the gapped value decreases the d.c. conductivity. Our calculations show that a decreases by 
increasing A as a function of n/rii. Moreover, the density dependence of a is linear at small d and A values and 
deviates from linearity at large d values/22 In the inset, we have shown the results for d = lOA which physically 
determine that the value of a increases by increasing d. Interestingly, a large value of a will be obtained for suspended 
graphene or by using the Si02 substrate instead of using boron nitride or SiC which result in the opening a gap due 
to symmetry breaking between sublattices. 

We also calculated the exchange and correlation energies as a function of A for various values of the cutoff A. 
The results are summarized in Fig. 4. we have found that the band gap effects become more appreciable at large 
cutoff values. The exchange energy is positive^ because our regularization procedure implicitly selects the chemical 
potential of undoped graphene as the half gap energy; doping either occupies quasiparticle states with energies larger 
than A, or empties quasiparticles with energies smaller than —A. Figure 4(b) shows the correlation energy fee as 
a function of A. Note that fcc has the same density dependence as 5ex apart from the weak dependence on A. In 
contrast to the exchange energy. Figure 4(a), the correlation energy is negativei^. It is important to note that there 
is a similar behavior between the kinetic energy and the exchange-correlation energy as a function of A. The kinetic 
energy is a slowly varying function to A at small gap values and increases by increasing A in middle and large values. 
Consequently, the total energy increases as a function of A. Fig. 5 is shown the total ground-state energy. In the 
inset, the total energy per particle is shown as a function of A for various values agr at A = 50. 

Figure 6 shows the charge compressibility, k/kq scaled by its noninteracting gapless compressibility as a function 
of A for different A values. The behavior of k suggests some novel physics qualitatively different from the physics 
known in the conventional 2D electron gasj ^^'^^ Kinetic energy and the exchange-correlation energy make negative 
contributions to the compressibility and therefore reduces the compressibility by increasing the A. 



IV. CONCLUSION 

We have studied the ground-state thermodynamic properties of a gapped graphcne sheet within the random phase 
approximation (RPA). Note that for a doped graphenc the Fermi hquid description is vahd. Our aim in this paper is 
investigating the ground-state properties of a gapped graphene sheet by going from a system with a hnear dispersion 
relation with vanishing the energy gap, A = to a system with a parabohc dispersion relation where A -^ oo. 
To achieve this goal, we have calculated the band gap dependence of noninteracting dynamic polarization function 
for doped graphcne sheet. As a consequence. We have presented results for the conductivity suppression over a 
wide range of energy gap. We have presented results of ground-state energies by incorporating many-body electron- 
electron interactions via RPA for gapped graphcne sheet. The total ground-sate energy increases by increasing the 
band gap values. This manner occurs based on our model Hamiltonian. We have finally presented results for the 
charge compressibility suppression over the energy gap. Importantly, the impact of gap energy on the thermodynamic 
properties would be noticeable for A > 0.2eF. 

Our results demonstrate the importance of including correlation effects together with the gap effects in the ther- 
modynamic quantities of a gapped graphcne. It should be possible to extend our work to include disorder effects. 
Another direction would be to consider the effects of temperature in the thermodynamic quantities. 
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FIG. 1: (Color online) Noninteracting dynamic polarization function for both gapped and gapless graphene in units of density 
of state, -D(eF) as functions of q/kp and uj/sf- 




FIG. 2: (Color online) (a); Static noninteracting polarization function as a function of g/kp for various A. (b): X isi 
fcp,itj,/x) as a function of lu/ef for various A. 
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FIG. 3: (Color online) (a): Conductivity as a function of njui for several energy gaps at d = 1 A, in the inset at d = 10 A(b): 
Conductivity as a function of A for various electron densities per impurity density. 
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FIG. 4: (Color online) Exchange (a) and correlation (b) energies as a function of A for various cutoff A at agr = 2. 
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FIG. 5: (Color online) Total ground-state energy per particle as a function of gap energy, A for various values of the cutoff 
A at a.gr = 2. In the inset the total ground-state energy per particle is shown as a function of A for various values of agr at 
A = 50. 
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FIG. 6: (Color online) Compressibility k/kq sacled by that of a noninteracting gapless system as a function of A for various 



agr and A. 



